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Abstract. The stellar initial mass function (IMF) determines the relative number of 
stars born at a given mass. Despite the tremendous effort to establish a universal IMF, 
the astronomical literature offers a wealth of diverse evidence for IMF variations. This 
review was prepared for a controversial debate at the conference "Starbursts - Near 
and Far" at Ringberg Castle, 2000, and gives a one-sided portrayal in favour of IMF 
variations. I will summarise the empirical evidence that the IMF varies with time, with 
environment, and for all stellar masses. While I see no obvious systematic trend in most 
regions of our Galaxy, there is at least an indication that the IMF is biased towards 
more massive stars in the early universe and in starbursts. 

1 Introduction 

The masses of stars span at least 3 orders of magnitude, from approximately 100 
solar masses for the highest mass stars to the hydrogen burning limit of 0.1 solar 
masses. The relative numbers of stars born at a given mass is described by the 
initial mass function (IMF). This distribution is extremely important for many 
fields in astronomy, from the theory of star formation to the interpretation of 
integrated properties of galaxies at the highest redshift. Although many inves- 
tigators have measured the IMF in a variety of star forming regions, we have 
not yet been able to come up with a final conclusion about the universality of 
the IMF. Specifically, Scalo's article |58| for the recent conference "The Stellar 
Initial Mass Function" has revived the discussion, and many scientists - 
including the author — have not decided on their final opinion. Nevertheless the 
conference organisers have chosen to split the review on the universality of the 
IMF to keep the discussion lively and controversial. This contribution focuses 
strictly on the empirical evidence in favour of IMF variations. I will not try to 
balance any arguments, but leave the critical comparison with the evidence in 
favour of a universal IMF || to the reader. 

The postulation of a universal IMF — in its strongest form — states that 
the IMF is and has always been the same in all regions of star formation every- 
where; the frequencies of initial stellar masses in any unbiased sample are always 
drawn from the same statistical distribution. I will present evidence against this 
postulation with respect to "is and has always been" as well as to "everywhere" , 
and I will also argue against this postulation being valid even in star forming 
regions of our own Galaxy. 

Although challenging the concept of a universal IMF, I refer to a reference 
IMF for comparison. This IMF (small inset in figure [j]) is the field star IMF 
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defined by Scalo |37j i n 1986, because its qualitative behaviour seems to be 
typical in many respects: (1) For high stellar masses the IMF can be described 
by a power law d f a ^ M oc M r . The exponent is referred to as slope of the 
IMF. In the notation of this article the value for the Salpeter |35| IMF is -1.35. 
This power law results from scale free star formation processes, for example 
dominated by turbulent pressure. (2) The IMF is flat — the slope of the power 
law approximately — for masses around 0.5 M©, indicating a characteristic 
stellar mass, for example the thermal Jeans mass. (3) The IMF declines for low 
stellar masses and may be described by a log-normal distribution, indicating a 
large number of independent parameters in the star formation process. However, 
I will present evidence that the IMF differs in various regions from this reference 
IMF in shape (e.g. log-normal versus power law and turnover versus continuous), 
in slope (e.g. bias towards high or low masses) and in characteristic mass. 



2 The Early Universe 

If the IMF depends on any environmental conditions, one would expect the 
largest deviations for the most extreme environments. The early Universe is 
such an environment. There are two very compelling indications that the IMF 
was different in the early Universe: First, no metal-free stars and only few very- 
low-metallicity stars have ever been found [fi]||], although low-mass stars, which 
have formed in the metal-poor early Universe, should have survived the last 
15 billion years. Only ss 500 stars are known with [Fe/H] « [Fc/H]q — the 
lowest metallicity observed in Galactic globular clusters — , and only w 100 stars 
with [Fe/H] w [Fe/H]©. No single star with primordial metallicity [Fe/H] 
w 100 oo 00 [Fe/H]© has ever been observed. Second, the solar neighbourhood is 
deficient in metal-poor stars ||(|. This discrepancy with the field star IMF is 
known as the "classical G-dwarf problem" . Because the lifetime of these stars is 
longer than the age of the Universe, we would expect to find many more such 
stars, if low-mass stars have formed with the same frequency in the early days 
of our Galaxy as they do now. As outlined by Larson [Q, this evidence finds 
its natural explanation in a time dependent IMF with a bias towards massive 
stars in the early Universe. Such a varying IMF is favoured for other reasons 
as well. For example, the gas of galaxy clusters is very hot and contains a large 
mass of heavy elements, which can not be explained from the proportionally 
few massive stars of the present day IMF. Also the observed strong evolution 
of the cosmic luminosity density with redshift is easier to explain with a top- 
heavy IMF. And microlensing experiments indicate that a significant fraction of 
dark matter is in the form of stellar remnants, which can only be produced from 
massive progenitors. See Larson's article JlJ] for a more detailed explanation 
and references. 
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Fig. 1. Evidence for IMF variations: This figure shows a comparison of the IMF of 
several selected Galactic regions. Variations of the IMF are evident for all stellar masses: 
a) The exponent of the IMF differs by more than one in the high mass star forming 
regions NGC 3603 |§ and Upper Scorpius OB2 §]. b) Compared to the Pleiades 
the IMF of the Orion nebular cluster |n| is deficit in stars with masses between 0.5 
- 1.25 Mq. c) The low-mass IMF in globular clusters is fitted best by a log-normal 
distribution but the field stars follow a power law in the same mass interval JS^]. 
The small inset shows the field star IMF of Scalo lB7| for comparison. 
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3 Starbursts and Extragalactic Star Forming Regions 
3.1 Starburst Galaxies 

Starburst galaxies are another extreme environment for star formation. Their 
star formation rate is several magnitudes larger than in our own Galaxy. So if 
the thermal Jeans mass is of any importance for a characteristic stellar mass, 
the heating from the unusually strong star formation should give preference to 
more massive stars. 

The starburst galaxy M 82 has been the prime example for a top heavy 
IMF for 20 years jn). As the closest starburst galaxy with a distance of only 3.2 
Mpc, no other starburst galaxy has been observed in such detail. The appearance 
of this galaxy is dominated by its nuclear starburst, which has probably been 
induced by a close encounter with M 81 about 100 Myr ago. The basic properties 
of the nuclear starburst are J24[|: Bolometric luminosity > 5.0 x 10 10 L©, absolute 
K magnitude < —22.5, ionising radiation Lt, yc > 10 54 photons/ s, total mass in 
the starburst < 2.5 x 10 s M Q , CO index > 0.21 and supernova rate w 0.1 



yr~ . Rieke et al. [pl|j32| have pioneered the detailed modelling of this starburst 
galaxy, and have come repeatedly to the conclusion that stars with masses below 
a few Mq form much less often in M 82 than in the solar neighbourhood. The 
argument for a top heavy IMF is basically that the infrared brightness — tracing 
the bolometric luminosity — is too high for the (dynamically measured) stellar 
mass, if one assumes a IMF similar to the field star IMF of Scalo p7|| . The 
ionising radiation, the CO index and the supernova rate provide the necessary 
constraints for the age of the stellar population. 

Although recent modelling fl36p[] with a more complex spatial structure and 
temporal evolution can not exclude a Salpeter power-law IMF extending down 
to subsolar masses, the starburst galaxy M 82 remains one of the most cited 
examples for IMF variations. 



3.2 Extragalactic Super Star Clusters 

The same technique can also be applied to extragalactic super star clusters. 
Spectral diagnostics, such as hydrogen emission lines and stellar absorption fea- 
tures, enable us to date the starburst. Dynamical mass estimates together with 
infrared observations provide a reliable mass-to-light ratio, and these together 
constrain the shape of the IMF. Compared with starburst galaxies, individual 
super star clusters offer the big advantage of a simple spatial (only one clus- 
ter) and temporal (only one burst) structure. Nevertheless these clusters can be 
looked upon as the building blocks of starburst galaxies, and should thus trace 
the IMF for rather extreme environmental conditions, too. However, the IMF in 
starburst clusters is not necessarily biased towards massive stars. 

The first super star clusters with dynamical mass estimates were NGC 1569A 
and NGC 1705-1 @. With a mass of 1.1 x 10 6 M (NGC 1569A) and 2.7 x 10 5 
Mq (NGC 1705-1), these clusters are both gravitationally bound and could 
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evolve into globular clusters p2[. However, recent evolutionary synthesis mod- 
elling of the mass-to-light ratio ]39| revealed that the IMF must be very different 
in the two regions. While the IMF in NGC 1596A is steep (slope ss -1.5) and 
extends to below 0.1 M©, the IMF in NGC 1705-1 is either flat (slope > -1) or 
truncated at a lower mass limit of 1 - 3 Mq. 

A similar conclusion has been found for the Antennae ]25} ]. The new gener- 
ation of infrared spectrometers has allowed the dynamical mass determination 
of young embedded super star clusters for the first time. The authors [^5| have 
compared each mass measurement with the results with the stellar synthesis 
model Starburst99 Q, that fit best the observed NIR- luminosity, Br7, and the 
CO and Ca absorption features. The measurements in two of the observed clus- 
ters are consistent with a Salpeter IMF extending to subsolar masses. However, 
in their cluster #2 [ft3| , such a power law IMF with a slope of -1.35 would in- 
clude all the dynamical mass (« 1.6 x 10 6 Mq) in stars more massive than 1 
Mq. Therefore no mass is left for subsolar mass stars. This cluster seem to have 
formed proportionally more massive stars. 

In conclusion, some of the measured mass-to- light ratios in super star clusters 
support the theoretical arguments for a top heavy IMF in starbursts, and thus 
systematic variations in the IMF, but the same kind of measurements also point 
to random variations in the IMF from cluster to cluster, even those in similar 
environments. 



3.3 Extragalactic Star Formation Complexes 

Unfortunately there are only a few dynamical mass measurements for star clus- 
ters. Statistics on IMF variations for a larger number of extragalactic star form- 
ing regions has to rely on integrated photometric and spectroscopic properties. 

A recent example for such a statistical approach is |$4[, who have included 
105 extragalactic star formation complexes for which UBVR photometry, Lyman 
continuum flux, metallicity and extinction measurements were available. These 
properties trace essentially the high mass stellar content, and the results are 
restricted to the IMF for stars more massive than 10 Mq. The models include the 
slope of the IMF, the maximum stellar mass and the age of the cluster. The star 
formation history is assumed either as a delta burst or as continuous. The average 
slope of the IMF was found to be -1.42. The standard deviation in the measured 
distribution is 0.91, and the authors estimate the standard error of their method 
to be 0.51. Therefore the measured distribution of slopes is meaningful f34| . 
Because the total number of analysed star formation complexes is rather large, 
the formal probability of measuring such a large standard deviation is basically 
negligible (< 10~ 10 ), if one assumes a universal IMF. 

Another investigation [ p3| has concentrated on 17 star formation regions in 
a sample of blue compact and irregular galaxies, using the integrated properties 
from UV, optical, FIR and radio measurements. This analysis uses the equivalent 
widths of H/3, Si IV, and Ca IV to constrain the relative number of 03-08 stars to 
B0-B3 stars, and thus the slope of the high mass IMF for stars more massive than 
10 Mq. The evolutionary tracks of starburst regions with three modelled IMF 
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slopes of 0.0, -1.35, -2.0 are clearly separated in a W(H/3) over W(Si IV) /W(C IV) 
diagram. The derived slope for the starburst regions varies between (3 regions), 
-1.35 (8 regions) and -2.0 (5 regions). Despite the variation in the derived slopes 
of the IMF, the authors interpret their result as in supporting a universal IMF, 
specifically because they see no trend with metallicity |2j|. However, the scatter 
in the diagram is much larger than the average error bars, and I would interpret 
the same diagram as evidence for a varying IMF. In addition, the authors also 
report that objects showing flatter IMFs are always small compact star forming 
regions, that four of the galaxies with optical continuum dominated by a previous 
generation of stars have IMF slopes close to -2, and that a previous burst of star 
formation could have hampered the formation of lower mass stars. 

4 IMF Variations in the Galaxy 

A direct measurement of the IMF from star counts is still restricted to our Galaxy 
and the most nearby galaxies. While we can only compare the relative number 
of high- to low-mass stars or constrain the slope of the high-mass IMF in the 
early Universe and in extragalactic starburst regions, we are able to trace the 
IMF down to substellar masses in nearby star forming regions. Figure |] shows 
a comparison of several IMFs measured in the Galaxy. We find variations in the 
slope and/or the shape of the IMF for all mass ranges. 

4.1 The Most Massive Star Forming Regions 

The most massive stars found in the Galaxy have about 100 M . Most of these 
stars have been found in either very compact star clusters or rather loose OB 
associations. The observations of these regions have improved significantly in 
the last few years with the HIPPARCOS mission, the Hubble Space Telescope 
(HST) and adaptive optics assisted telescopes. 

For example, HIPPARCOS has revealed 178 new members of the Upper Scor- 
pius OB2 association Only 91 stars were previously known in this association. 
Figure [j] shows the preliminary IMF. With a slope of w -1.9, this IMF resembles 
very well Scalo's |57j field star IMF. I have compared this high-mass IMF with 
the IMF in NGC 3603, the most massive visible HII region in our Galaxy 
Adaptive optics assisted observations have revealed more than 800 stars in the 
central parsec of this cluster. Although the IMF shows no turnover or truncation 
down to < 1 Mq, its slope is only « -0.7 for stars with masses between 3 and 30 
Mq. A very similar result was found with the help of HST in the Arches cluster 
near the Galactic centre ||. The slope of the IMF for stars more massive than 
10 M Q is w -0.65, much shallower than the reference field star IMF, too. The 
Arches cluster and NGC 3603 are two of the most massive young clusters in the 
Galaxy, and each cluster has a few 1000 Mq in O-stars. Therefore the difference 
between the field star IMF and the IMF in NGC 3603 and the Arches cluster 
can not be explained by small number statistics. Such an argument is often used 
to invalidate the evidence for IMF variations from star counts |D| . The finding 
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in these two clusters supports the hypothesis that the IMF in starburst regions 
is biased towards massive stars. However, the deficit of low-mass stars is caused 
by a shallow IMF, and does not result from a truncated IMF, which has been 
indicated in M 82 

4.2 Associations and Open Clusters 

The discussion about a universal IMF was revived by Scalo p£| in 1997, at a 
time, when the majority of the participants of a conference on the stellar initial 
mass function was arguing for a universal IMF He made his case for IMF 
variations basically from a statistical interpretation of the measured slopes of 
the IMF in Galactic and Magellanic cloud associations and open clusters. Scalo 
plotted the slope of the IMF of 61 clusters against the average logarithmic stellar 
mass considered in the measurements. Although the typical error in the slopes 
is 0.1 - 0.4, the spread in the diagram above w 1 M is so large that the author 
saw no basis for adopting some average value. 

A serious argument against such evidence for IMF variations is that the 
results for different clusters have been obtained from different authors using dif- 
ferent techniques. Therefore the IMF variations may have been mimicked by 
systematic errors, which have not been addressed carefully enough. I have thus 
repeated a similar x 2 analysis for this conference to estimate how large the sys- 
tematic errors would have to be, if the IMF were to be universal. My compilation 
includes the slopes of the IMFs for 51 clusters with intermediate- to high-mass 
star formation from 7 publications. The publications have been selected such that 
each paper includes more than 4 clusters, so that I could carry out statistics for 
each author and technique separately. Table |l| summerises the IMF properties 
of the different publications. Clusters with suspicious IMF measurements (for 
example NGC 436 in [^], which has probably undergone dynamical evolution), 
and results without error estimates (for example NGC 7235 in [|l)) or with a dif- 
ferent mass range (for example NGC 376 and N 24 in [g) have not been included 
in the statistics. For the ease of interpretation, I have converted the reduced % 2 
values into a significance level |HJ , which is the probability of measuring such a 
high or higher x 2 , if the IMF would have a universal slope equal to the weighted 
mean slope in this subsample. Even at a significance level of only 5 %, we have 
to reject the hypothesis of a universal IMF in 5 of 7 cases. The statistics for the 
whole sample is even worse. The standard deviation in the slope distribution is 
1.9 times the mean error in the measurements of the slopes. The likelihood for 
such a large discrepancy is basically negligible (< 10~ 10 ) for a universal IMF. If 
the IMF is indeed universal, the authors must have underestimated their errors 
by a factor of two. The systematic errors would have to be even larger if we 
include the field stars of the Magellanic clouds, which show a particularly steep 
IMF with a slope of w -3 to -4 |f|. 

A very intriguing example with very different IMFs at intermediate stellar 
masses are the Pleiades and the Orion nebular cluster. These clusters have been 
the target of some of the most careful recent measurements of the IMF |2(].|l(J . 
The IMFs of these two clusters are also displayed in figure |[ In contrast to the 
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Table 1. IMF variations in 57 clusters with intermediate- to high- mass star formation 



Ref. Number 


Weighted 


Standard 


Mean error 


Reduced 


Significance 




of 


mean 


deviation 


in slope 


x 2 






clusters 


slope 


in slopes 


measurements 








40 


4 


1.43 


0.34 


0.15 


5.15 


0.15 % 




33 


5 


1.27 


0.80 


0.33 


14.82 


« a 




27 


6 


1.22 


0.31 


0.30 


2.92 


1.2 % 




21 


12 


1.12 


0.44 


0.34 


1.78 


5.3 % 




22 


5 


1.27 


0.18 


0.16 


1.64 


16.2 % 






12 


1.92 


0.54 


0.34 


3.22 


0.021 % 




29 


7 


1.42 


0.26 


0.15 


18.08 


« a 


all 


51 


1.52 


0.53 


0.28 


6.35 


fa a 



a The significance level is < 10 



IMF of the Pleiades, which follows the field star IMF, the IMF in Orion shows a 
regime of "missing" stellar masses between 0.5 - 1.25 Mq. Other examples are 
NGC 6231 and NGC 2264. The IMFs in these clusters have been measured in 
both cases with the same technique — UB VRI and Ha photometry — and by the 



same authors |41,42|, therefore minimising systematic variations. Nevertheless, 
the IMF in NGC 2264 rises continously down to below 0.6 Mq, whereas the IMF 
decreases abruptly below 2.5 M Q in NGC 6231. 



4.3 Low-Mass Stars in Young Clusters, Globular Clusters and the 
Field 

If we assume a Salpeter (35) power-law IMF with a slope of -1.35 extending down 
to a lower mass limit of 0.1 Mq, 96 % of all stars would have a mass smaller 
than 1 Mq , and 55 % of the stellar mass would be included in stars less massive 
than the sun. However, these low-mass stars are very faint, and many such stars 
have not been detected before the large ground based infrared surveys and HST. 
These observations have also revealed IMF variations in many regions. 

To illustrate these variations, I will again compare different regions with the 
field star population. Figure [l] shows a recent measurement of the field star 
IMF for low-mass stars. It was derived from the Deep Near-Infrared Survey (DE- 
NIS) and the 2 Micron All-Sky Survey (2MASS). In the 0.1 - 1 M mass range, 
this IMF can be represented by a power-law mass function with a slope of « 0.1. 
The statistical uncertainties are « 0.13. A similar IMF has been observed in the 
young star forming regions IC 348, p Ophiuchi, and the Trapezium Jl6|Jl~^JTc| . 
Their IMF is flat or slowly rising — slope J> — from the brown dwarf regime 
to 0.6 — 1 Mq, where it rolls over to a power-law with a slope of w —1.7. It is 
important to note that the IMF is not log-normal |Q . 

In contrast, recent HST observations indicate that the IMF of globular clus- 
ters can be described by a log-normal function with a peak near 0.3 Mq pq. 
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Figure [j] shows the best fit for a sample of 12 globular clusters. The authors p8[ 
explicitly exclude a single power-law IMF in the 0.1 — 0.6 Mq range. Therefore 
the IMF of globular clusters seems to be fundamentally different from the IMF 
in young star forming regions and the field. 

A meaningful interpretation of the IMF measurements for masses below the 
hydrogen burning limit is increasingly more difficult, because the available data 
are sparse. However, observations of young star forming regions indicate IMF 
variations even in the mass regime of brown dwarfs. The Taurus star forming 
region, for example, is significantly deficient in objects below 0.1 Mq compared 
to the Trapezium fl9| . If the IMF of both objects is normalised by the number 
of stars between 0.1 - 1 Mq, then w 13 brown dwarfs with masses > 0.02 M 
are found in the Trapezium, but only one in Taurus. 

5 Conclusions 

Variations of the stellar initial mass function have been reported for all masses 
and in a large variety of stellar populations. We find evidence for IMF varia- 
tions in the early Universe, in starburst galaxies and extragalactic star forming 
regions, in Galactic star clusters and associations, and in the field star popu- 
lation. Neither the shape nor the slope and the characteristic mass of the IMF 
seem to be excluded from these variations. There is some indication that the 
IMF is systematically biased towards more massive stars in the early Universe 
and in starbursts. However, I see no obvious systematic trend in those regions 
where the IMF could be constructed from direct star counts. If the IMF is indeed 
universal, it will be very difficult to prove this postulation empirically, because 
every individual measurement of a different IMF has to be invalidated. Only 
a consistent theory of star formation, with clear and testable predictions, will 
finally convince in the scientific world. 
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